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AN ANALYSIS OF THE OPERATION OF A MASS MEASURING SYSTEM 
IN AN ORBITING SPACECRAFT 

v By Bruce A. Conway 

Langley Research Center 

SUMMARY 

An analysis of the operation of a mass measuring system, which is based on the 
oscillating spring -mass principle, aboard a spacecraft is made. This analysis considers 
the dynamic effects of the oscillating device on an uncontrolled spacecraft by developing 
the equations of motion for the system. An approximate solution of these equations is 
then deduced, and a simple analytic expression for the maximum attitude error is found. 

Numerical integration of the equations of motion was performed on a digital com- 
puter for some typical measurement operations. Comparison between the approximate 
and exact solutions for these cases was very close. Results showed small (less than 
0.02°) attitude errors arising from the operation of this device. 

INTRODUCTION 

With the advent of long-term spaceflight (from 2 weeks to a year or more) it 
becomes necessary to monitor various physiological indicators of an astronaut's well- 
being. One important factor is his weight, which can serve as a monitor for other func- 
tions (such as cardiovascular effort or muscle condition). It then is important to find 
some means of determining astronaut weight in the zero-gravity environment of manned 
spacecraft. Studies have shown that an oscillating spring-mass system is perhaps an 
optimum method of determining the unknown mass of astronauts and other spacecraft 
apparatus (ref. 1). Life support expendables (e.g., food) and waste products (CO2 cani- 
sters and urine samples) could be weighed by using this technique. The required spring- 
iness system could be lightweight and easily operable by a single crew member. 

A prototype mass measuring system (MMS), employing the oscillating spring-mass 
principle, has been designed and fabricated under NASA contract. (See ref. 2 and fig. 1.) 
This report discusses the effects of operating a device similar to the prototype aboard an 
uncontrolled spacecraft with results presented for some typical operations. 


SYMBOLS 


C^Cg constants of integration 

F X’ F Y’ F Z force along X-, Y-, and Z-axis, respectively 

Hem angular momentum of spacecraft about its center of mass 

H 0 angular momentum of spacecraft in geometric coordinates 

H s angular momentum of spacecraft mass center 

[i] inertia tensor 

^X’ly^Z spacecraft moments of inertia 

^XY’^XZ’^YZ spacecraft products of inertia 
K spring constant 

M X’ M Y’ M Z moment about X-, Y-, and Z-axis, respectively 
m single moving mass 

mj arbitrary mass 

m g spacecraft mass 

Q mass factor 

Rq radius of geometric origin with respect to inertial axes 

r radius to mass measuring system 

r mj’Hmj radii of arbitrary mass 

r ms’Hms radii of spacecraft mass center 

Tx,Ty,T2 torque about X-, Y-, and Z-axis, respectively 
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t 


time 


X b> Y b> Z b body-fixed spacecraft axes 

Xj,Yj,Zj inertial axes 

moving axes parallel to inertial set 
x,y,z rectangular coordinates 

ot X9 oty,a z attitude angles of spacecraft obtained from approximate solution 

Ar vector representing oscillation of mass measuring system 

Ar x ,Ary,Ar z amplitude of mass measuring system oscillation in x-, y-, and 
z -direction, respectively 

(p 9 6, \f/ modified Euler angles 

tt> n natural frequency of spring-mass system 

oo X9 <*)y 9 u) z angular velocities of spacecraft about body axis 

Subscripts: 

c quantity resulting from cocking of mass measuring system 

o referenced to geometric axes 

s referenced to spacecraft mass center 

Dots over symbols indicate derivatives with respect to time. 

A bar over a symbol indicates a vector quantity. 

ANALYSIS 

Measurement Technique 

The principle of the oscillating spring-mass system is applied by placing the 
unknown subject mass in a carriage that is constrained (mounted on a rail, for example) 
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to move in one degree of freedom. The carriage is connected to the spacecraft structure 
by two springs in parallel (one at each end of the carriage). The carriage, with its sub- 
ject mass, is displaced from its equilibrium position (cocked) and then released to oscil- 
late at the natural frequency of the carriage- spacecraft system. By measuring the fre- 
quency (or the period) of oscillation and by knowing the spring constant of the springs the 
unknown mass can be determined. Consider the following derivation. 

From figure 2, the following expressions are evident: 

mjXj + Kx = 0 (1) 

m 2 x 2 - Kx = 0 (2) 

where spring displacement is 

x = Xj - X£ (3) 

Multiplying equation (1) by m 2 and equation (2) by mj and subtracting equa- 
tion (2) from equation (1) give 

m^rngX + K(ni2 + m-^x = 0 (4) 


m 0 + m 1 

x + K — - lx = 0 

m l m 2 


(5) 


vibrations having a natural fre- 

( 6 ) 

For the present system (see appendix) 

m 1 = m (7) 

m 2 = m s - m (8) 


Equation (5) is the equation of motion for harmonic 
quency <u n of 


w n 




[ K^m 2 + rrtjj 


m l m 2 


and 


m l m 2 
m 2 + m l 


m(m s - m) 


= Q 


(9) 


Equation (6) can thus be written 



( 10 ) 
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The symbol w n is the frequency that would be measured if there were no non- 
linearities in the system. However, friction will always be present to a certain extent, 
and spacecraft body rates will probably exist (introducing additional degrees of freedom). 
In the following analyses, these rates are assumed small, so that a one-degree-of-freedom 
system is maintained. In addition friction is assumed negligible. In actual space-flight 
use, empirical calibration would be used to reduce measurement data. 


Applied Torques 


The disturbance moments M 0 (see appendix) resulting from operation of the mass 
measuring system are caused by the forces exerted by the springs. These spring forces 
are 


F X = K(x - x 0 jj 
F y = K(y - y Q ) ) 
F z = K(z - z 0 ) 


(ID 


where x 0 ,y 0 ,z 0 are body coordinates for the equilibrium position of the device. The 
x, y, and z coordinates for the moving mass are given by (assuming that velocity of 
the MMS following release is in a positive direction) 


x = x Q - Ar x cos o> n t 
y = y 0 - Ar y cos w n t\ 


( 12 ) 


z = z 0 - Ar z cos o) n tj 


The components of M 0 due to the spring forces are written 


M x = yFz " zF Y 
My = zF x - xF z > 


(13) 


M z = xFy - yFxj 

By using equations (11) and (12), equations (13) may be rewritten as 

M x = K cos w n t^z 0 Ary - 2Ar y Ar z cos o> n t - y 0 Ar z ^ 

My = K cos o> n t^x 0 Ar z - 2Ar z Ar x cos u> n t - z 0 Ar x n 

M z = K cos u> n t(y 0 Ar x - 2Ar x Ar y cos o> n t - x Q Ar y ) 


(14) 


Substitution of equations (14) into the set (A39) yields the spacecraft equations of motion 
for the present problem. Other disturbances such as gravity gradients and aerodynamic 
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torques have been neglected since they cause negligible rate and attitude changes in the 
short time (10 to 15 sec) needed for a measurement operation. 


Approximate Solution 


A linearized solution to the equations of motion can be obtained, which will prove 
useful in many applications not requiring the exact solution. If it is assumed that pro- 
ducts of inertia are negligible (=0), then the equations of motion are Euler's equations for 
principal axes, that is 

p]aJ+ cox [l]aJ = T (15) 


or, expanding and arranging terms 


T X % - 

co x = -y— + — — w y w z 


J X : X 


• Ty ^7 - ly 
Wy = — + — ^ W X U> Z 


!y Iy 


u) z = — + 


T Z *X " l Y 

i z + ”i^“ WxW l 


(16) 


If small angles and angular rates are assumed ( a <0.1 rad, cu < 0.1 rad/sec) then the 
terms containing products of angular velocities will be negligible in comparison with the 
T/l terms and the equations of motion can be written as 



[I] w = T 

(17) 

Accordingly, rewriting equation (17) for the mass measurement application gives 



[iJcj + ? X Qr = M 

(18) 

But 


M = ? X F 

(19) 

and 


F = K(r - r 0 ) 

(20) 

so that 


= p]’ 1 ? x K(? - ? 0 ) - QrJ 

(21) 

Now 


r = r 0 + (cos 0 ) n tjAr 

(22) 


r = -fcu n ^cos o) n t^Ar 

(23) 
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Substituting equations (22) and (23) into equation (21) and collecting terms gives 

co = {TJ -1 r 0 x |(k + Qco n 2 ) cos co n tjAr 
Integrating twice to find rate co and attitude a results in 


( 24 ) 


w = c i r lf o x 


+ Qco n 2 


CO 


n 


sin co n yAr + Cj 


(25) 


and (since angular rotations can be treated as vectors for small angles) 


a = -[lj _1 f 0 


(K + Qco n 2 


CO 


cos co n t Ar + Cjt + C 2 


(26) 


n 


At t = 0 for any measurement operation co = 0 and a - a a (static angle caused by 
cocking the MMS), and the constants of integration and C 2 can be evaluated as 

follows: 


c x = 0 


c 2 = pcr lf ox 


K + Qco n 2 


Ar + a. 


co 


n 


Equations (25) and (26) can now be written as 

. (k + Qco 2 


CO 


n 


a = -H’^o 


X 


K + Qco 2 


co r 


ntjAr 

(27) 

,1_ _ 

lj Ar + o-o 

(28) 


In cocking the MMS, the cocking force is equal to or greater than the spring force, so that 

co c = (X)” 1 ? x (-Qr c ) (29) 


where 


r = r n + r- 
o r 


a o = -QP]" 1 Jj( ? o x f c + x ? cj dt 


(30) 


The cocking acceleration vector r c is collinear with the incremental cocking radius 


r-* and hence 
r 
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r u x r c = 0 
r 


a o = -QP]" 1 jj ? o x f c dt = -Q[I] _1 ? 0 x Ar 
Equation (28) can now be written to include equation (31) 

“--prV 


K + Qo) 2 

;^-(cos “n 1 ” X ) + Q 


Ar 


(31) 


(32) 


Orientation of MMS to Yield Minimum Disturbances 


Since small spacecraft attitude angles have been assumed, a can be treated as a 
vector, and solutions for oj and a can be expressed in terms of the components of 

03. f o, and Ar, as 


co = 


K + Qco n ^ 


w n 

y ( 


IlT 


+ ^( r O. 


Ar,, - r Ar_V 
x y o,y x ) 


1/2 


sin co n t 


(33) 



For any time t the magnitude of a and co for a measurement operation is dependent 
upon the quantity 



o,y 


Ar z 


o,z Ar y) 2 + A( r o,z ir x - r 0 ,^zf * ~ 2 ( r o,x Ar y 



2 


in equations (33) and (34). This quantity is represented by N in subsequent equations. 
This quantity is a function of the spacecraft moments of inertia and the location and 
orientation of the MMS within the spacecraft. For a given spacecraft the moments of 
inertia are constant, and rate and attitude errors are functions only of position of the 
MMS and the orientation of its line of action. 
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If operation of the mass measurement device is assumed in a plane perpendicular 
to one of the spacecraft axes, then one component of Ar is zero. Accordingly, let the 
plane of operation be perpendicular to the spacecraft X-axis, so that Ar can be 
expressed as 


Ar = < Ar y > 


||at 


2 



(35) 


It can be seen from equations (33) and (34) that o> and a are minimized whenever the 
quantity N is a minimum. With the use of equations (35), 


N = 


ix 2 


r o,y 2 ( Ar2 ' Ar y 2 ) - 2 r o,y r o,z Ar y ( Ar2 “ Ar y 2 ) 1/2 + r o,z 2 Ar y 2 


+ 



o,x 


Ar z 






J Ar, ; 


(36) 


If a specific location within the spacecraft is considered, then r Q x , r o,y> anc * r o,z 
are constant, and the minimum values of a o and a are influenced by the magnitude of 
Ary (since Ar has a constant magnitude) or by the orientation of the mass measure- 
ment device with respect to the y and z axes (since Ar x = 0). This orientation can 
be specified by the angle 77 whose cosine is Ar^/Ar or 


77 = cos 


-1 



(37) 


To find the value of Ar y (and, hence, the value of 77 ) which will minimize N, take 


9N 

9Ar y 


= 0 


(38) 


and solve for A r^/A r . To obtain a closed-form solution, the binomial expansion of 

^Ar^ - Ary ^^ 2 must be used; this results in a ( 2 n)th degree algebraic equation (where 
n is the number of terms used in the expansion). A digital computer was used to solve 
this equation for cos 77 as a parametric function of r o x , io,y> and r 0?z . Figure 3 
presents cos 77 as a function of MMS location within the spacecraft for a typical inertia 
distribution g - h = 0.1 where 


g 




(39) 
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Digital Computer Simulation and Results for Typical Operations 

A digital computer program has been written that includes equations (A2) and 
(A39) along with equations (11), (12), and (13). This program employs a fourth-order 
Runge-Kutta integration scheme to integrate equations (A39) and (A2) to yield spacecraft 
angular velocities and Euler angles, respectively. All initial rates and Euler angles, as 
well as inertia properties and MMS properties can be read into the program as initial 
conditions and changed at specified times throughout the program. Program output con- 
sists of the Euler angles cp, 6, and \f/, as functions of time, and o> x , o> y , and co Z) 
also as functions of time. In simulating the MMS operations to be described, all equa- 
tions were solved exactly. Initial spacecraft rates and Euler angles were set equal to 
zero (prior to cocking of the device). First examine the response (rigid body) of a space- 
craft during a typical mass measuring system operation. An Apollo Applications 
Program spacecraft, consisting of an Apollo command and service module docked with 
an S-IVB Orbital Workship/ Multiple Docking Adapter combination, is considered. This 
configuration, with the orientation of its body axes, is shown in figure 4. The MMS is 
assumed operating in the Orbital Workshop. Moments of inertia of the spacecraft about 
its center of mass are 

J X,o = 115 000 slug-ft 2 = 155 917 kg-m 2 
I Y o = 1 970 000 slug-ft 2 = 2 670 926 kg-m 2 
I z Q = 1 930 000 slug-ft 2 = 2 616 694 kg-m 2 

All products of inertia are taken to be zero. The coordinates of the MMS equilib- 
rium position are 

r QjX = 15.0 ft = 4.5720 m 
r Q?y = 2.0 ft = 0.6096 m 
r 0?z = 3.0 ft = 0.9144 m 

and orientation of the device is such that (with Ar = 0.5 ft = 0.1524 m) 

Ar x = 0 

Ar y = -0.3535 ft = -0.1077 m 
Ar z = 0.3535 ft = 0.1077 m 

Characteristics of the device are those of the final prototype design (ref. 2). The 
weighing of a 160-lb (72.544 kg) man was simulated (including cocking and carriage 
return) with results shown in figure 5. This figure gives the time histories of the Euler 
angles <p, 6, and i//. Also shown in this figure are the approximate solutions to the 
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equations of motion for this case - a x (~cp), oty(~6), and a z (~\p) - given by the dashed 
curves on the plots. It can be seen, from comparison of the exact and approximate solu- 
tions, that the assumptions made in the derivation of the approximate solution (small 
angular velocities, etc.) are valid; this solution provides an acceptable means of esti- 
mating spacecraft error. 

The operation just considered did i.ot take into account the optimum orientation of 
the mass measuring device within the spacecraft. For the given location (with g - h 
essentially zero) figure 3 can be used to determine the proper value of the orientation 
parameter cos 77 which in this case is 0.554. Thus, since Ar x = 0, Ary and Ar z 
can be found from the following equations: 

Ary = Ar cos 77 = 0.227 ft = 0.0692 m 

Ar z = Ar)/l - cos ^77 = 0.416 ft = 0.1268 m 

With this orientation, the spacecraft response is given in figure 6 for a similar operation 
as the previous case. Again, the approximate solution is shown to be in good agreement 
with the exact solution. It can be seen from comparison of figures 5 and 6 that 
reorienting the device significantly reduces spacecraft attitude excursions. 

CONCLUDING REMARKS 

This report summarizes the effects of operating a spring-mass device onboard a 
spacecraft for the purpose of determining mass. It should be recognized that this analy- 
sis considers only a rigid subject mass. From the results of typical man-weighing 
operations the frequency of oscillation was approximately 25 cycles per minute, which is 
sufficiently low so that a properly restrained subject will behave very nearly as a rigid 
body. 

Since the subject mass and spacecraft properties were arbitrary throughout the 
theoretical development, it is evident that the analysis presented herein can be applied 
to masses of any magnitude in any spacecraft configuration. Because mass measuring 
requirements aboard a manned spacecraft may range from a fraction of a gram (as in 
chemical analysis) to several hundred pounds (recoverable satellites, for example) the 
technique and theory presented in this report may be applied to a wide variety of mass 
measurement problems. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 4, 1968, 

127-53-16-05-23. 
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APPENDIX 


DERIVATION OF SPACECRAFT EQUATIONS OF MOTION 


The motion of the spacecraft is defined with respect to a set of moving body-fixed 
axes Xk,Yb,Zk and these body axes will then be related to space-fixed axes X^Y^Zj 
by means of a set of modified Euler angles (ref. 3). These modified Euler angles, shown 
in figure 7, result from three consecutive rotations. The first, about the Zj-axis, carries 
the Xj- and Yj-axes through an angle \}/ measured in a horizontal plane. The second 
rotation, about the new Yj-axis, carries the Xj- and Zj-axes through an angle 9, mea- 
sured in a vertical plane. The final rotation, about the new Xj-axis, takes the Yj- and 
Zj-axes through an angle cp, measured in an inclined plane, to give the X^-, Y^-, and 
Z b -axes. 

The modified Euler angles may be expressed mathematically in terms of the angular 
rates about the body axes. These rates, as a function of Euler angle derivatives, are 


o> x = cp - \p sin 9 

u>y = 6 cos cp + xjs sin cp cos 9 > 

u> z = \jy cos 9 cos cp - 9 sin cp j 


(Al) 


Solving for the Euler rotations gives 

cp = o> x + o>y sin (p tan 9 + a> z cos cp tan 9 s 
9 = o>y cos cp - co z sin cp ) 

i//= co z sec 9 cos cp + coy sec 9 sin cp ^ 


(A2) 


Integration of the set (A2) yields the Euler angles for the spacecraft. 

Before developing the equations of motion, the coordinate systems must be defined. 
Figure 8 illustrates the notations to be used. The Xj-, Yj-, and Zj-axes are fixed in 
inertial space. The Xj-, Yj-, and Zj-axes have their origin at the spacecraft while 
remaining parallel to the Xj-, Yj-, and Zj-axes. The Xb,Yk,Zk system is a set of body- 
fixed spacecraft axes to which all coordinates within the spacecraft are referenced. 

In order to obtain the Euler rates, then, it is first necessary to determine the space- 
craft body rates. In the following analysis, these body rates will be determined from 
momentum considerations (ref. 4). 
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The angular momentum of the spacecraft, in its geometric coordinates is defined 


as 


^ ^mj x m j R mj 


(A3) 


and the momentum of its mass center is 

H s = ^ms x m s^ms 

The total angular momentum about the mass center is then 


(A4) 


H cm = H 0 - H s (A5) 

or 

H C m = J ? mj x m j( R o + ^mj + « x f mj) - r ms X m s (R Q + r ms + SJ X ? ms ) (A6) 

Now, 


. _ I m i f mj _ I m j f mj 

r ms v— rn 

2>i 

Rewriting equation (A6) by using equation (A7) results in 

V- - ^ v . , 

H C m = ^ r mj m j X R 0 ^ X m s R 0 + 2 r mj x m j r mj ' r ms x m s r ms 


+ J f mj x m j( w x ? mj) - f ms X m s( w x %s) 


(A7) 


(A8) 


or 

H cm = ^ x mj f mj - r ms x m s ? ms + ^ r mj x mj(w x r mj ) - ? ms x m s (w x r ms ) 

(A9) 

For a rigid body, ? m j = 0 and r ms = 0, and the angular momentum for the rigid body 
H rb becomes 

H r b = J f mj x m jp x f mj) " f ms x m s(^ x ? ms ) (A10) 
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or, from elementary vector analysis 

H r b m j[( f mj ' ? mj)^ " ( f mj ' “) f mj] - m s[( f ms ' r ms )w - ( r ms ■ w)r ms ] 
Expanding equation (All) and writing in matrix form gives 



H rb,x 


" ix 

-I XY 

_I XZ 



H rb = 

H rb,y 

= 

_I XY 

% 

-I YZ 




^rb,z_ 


_ _I XZ 

-I YZ 

*Z 


0> z 


(All) 


(A12) 


Here, the moments of inertia are about axes passing through the center of mass of the 
spacecraft (by the parallel axis theorem). 

Again consider the origin of the geometric coordinates at the center of mass of the 
spacecraft without any moving masses, and the moment of inertia about the X centroidal 
axis becomes 

x x = j m j(yj 2 + z j 2 ) + J m i(y 3 2 + z j 2 ) - m s(y s 2 + z s 2 ) ( A13 ) 

J =1 v ^ ) K ^ J 

Spacecraft Inertia of Inertia shift to 

inertia without (n - K) moving mass center 

moving masses masses 

= Constant ! X,o 


Now, if it is assumed that there is only one moving mass m located at a distance r 
from the geometric origin, the moment of inertia about the X-axis (centroidal) can be 
written 

J X = Ix,o + m (y 2 + z2 ) " m s(y s 2 + z s 2 ) (A14) 


Since r g is now 


mr 
r s = ^ 


(A15) 


and 



x x = : x,o + Q(y 2 + z2 ) 


(A16) 

(A17) 
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where 


Q = 


m (m s - m) 
m s 


(A18) 


All inertias can be written in the same manner as follows: 

l x = ix,o + Q(y 2 + z 2 ) 

X Y = r Y,o + Q( x2 + z 2 ) 

: z = l z, o + q ( x2 + y 2 ) l 

*XY = ^,0 + 

ixz = ixz,o + Q( xz ) 

J YZ = !yZ,o + Qfr 2 ) ^ 


(A19) 


The total moment about the geometric origin is 

M 0 = ^ r m j X mjR m j (A20) 

The force equation, through the mass center, is 

F s = m s f ms (A21) 

The symbol F s represents the total external force vector acting through the mass 
center. From equation (A3) 


= Jl (Z rm i x (A22) 

"dt~ = dt 3 x m i Rm i + ^ r mj x m i Rm i ( A23 ) 

= I x m i(^o + %*) + J f m j x m i Rm i < A24 ) 

IT = Z m i ^ X *0 + 2 * mj x mjR mj (A25) 

dH o dr „ v- _ — 

m g — ^ X Rq + 2 ^ r m j X mjR m j (A26) 
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Similarly, from equation (A4) 


Also, from equation (A5) 


dn d r • *• 

= m s -^ x *0 + f ms x m s^ms 


dH cm dH 0 dH s 


dt dt dt 


dH 


cm 


dt 


= X f mj x m j R mj - r ms x m s R ms 


(A27) 


(A28) 


(A29) 


Now, by substituting for R ms from equation (A21), equation (A29) will become 

dlL 


‘cm 


dt 


=1 


r mj x m j R mj ” r ms x R s 


(A30) 


Referring back to equation (A20) and considering equation (A30) then 

dH r 


M 0 ^ + ^ms x R s 


(A3 1) 


The applied moment about the geometric axes will equal the derivative of the angular 
momentum about the mass center if no external forces act through the center of mass. 
For most rotational problems the external force contributions due to mass center shifts 
can be neglected. In this case 

dH r 


M n = 


‘cm 

dt 


(A3 2) 


Referring back to equation (A9) and remembering that by letting ? m j and ? ms both 
equal 0, an expression for H rb was obtained, and the angular momentum about the cen- 
ter of mass can be written as 


R cm = ^ f mj x m j f mj " f ms x m s?ms + H rb 


Again, only one moving mass m is considered and equation (A32) can be written 

dH r 


Mr, = 


‘cm 


dt 


dt V 


r x mr - r ms x m c r_ 0 + H 


s 1 ms + n rb) 


(A3 3) 


(A34) 


Substituting for r ms and r ms by using equation (A15) gives 

*rb 


Mo = 


rxmr-~x m r) + H T 


m 


s 


\ 


(A3 5) 
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“o '! 1 


m(m s - m) 

(f X r ) + H K 

m 

' rb 


M ° = a 


(A3 6) 
(A3 7) 

(A3 8) 


e [q(? X r)+ H rb ] 

Expanding equation (A37) results in the following equation: 

Mq = qJ? X r + co X (r X $ + H rb + CJ X Hrb 

Now, by substituting for H r b (and H rb ) by making use of equation (A22) the following 
components of M 0 about the X^-, Y^-, and Z^-axes through the geometric origin are 
obtained: 

M x = Q jjyz - zy) + [w y (xy - yx) - co z (zx - xz)jj 

+ (^x + w x^X " I XY c ^y " w y*XY ” ^xz^z ~ ^z^xz) 

+ |coy^I z w z " Ixz^x " ^Z^y) - a) z( I Y w y " l YZ w z ~ I XY w x)j 

= Q j^(zx - xz) + jw z (yz - zy) - a> x (xy - yx)jj 
+ ^Iy^y + Wyiy - Iyz^z - w z iy Z - IxY^x - ^x^Xy) } (A3 9) 

+ j^z^X^x " I XY co y “ r XZ w z) - w x(lz w z " IxZ w x “ ^Z^yJ] 

= Q f(xy - yx) + [u> x (zx - xz) - w y (yz - zy)jj 
+ ^I z d) z + co z i z - Ixz^x “ ^x^XZ “ Jyz^ “ ^y^Yz) 

+ ^(^y “ l YZ^z ~ I XY w x) " co y( I X co x “ IxY^y " ^Z^zj] 


Mt 


Mr 


If the moments MxMy,M z and coordinates for the moving mass x,y,z are expressed 
as functions of time, then equations (A39) can be integrated to find u> x , coy, and co z with 
numerical methods. The spacecraft rates, then, are used in equations (A2) to find the 
Euler angle rates. 
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Figure 3.- MMS orientation parameter cos n as a function of location within the spacecraft. 
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Figure 6.- Response of spacecraft to typical MMS operation (with optimum orientation). 
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